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Abstract 

This paper considers an algorithm for synthesis of opti- 
mal controllers for full information feedback. The syn- 
thesis procedure reduces to a single linear matrix inequal- 
ity which may be solved via established convex optimiza- 
tion algorithms. The computational cost of the opti- 
mization is investigated. It is demonstrated the problem 
dimension and corresponding matrices can become large 
for practical engineering problems. This algorithm repre- 
sents a process that is impractical for standard worksta- 
tions for large order systems. A flexible structure is pre- 
sented as a design example. Control synthesis requires 
several days on a workstation but may be solved in a 
reasonable amount of time using a Cray supercomputer. 

Introduction 

Flexible structures present interesting and difficult prob- 
lems to control engineers. Lightweight structural ele- 
ments give rise to lightly damped, often closely ^spaced 
natural frequencies. Linear, time-invariant mathemati- 
cal models of these structures do not accurately repre- 
sent the true system. Uncertainty is included to rep- 
resent unmodeled dynamics and high frequency modes. 
Parametric uncertainty may also be included to repre- 
sent unknown coefficients in the state-space models such 
as natural frequency and damping coefficients. 

There are many control design techniques for synthesis of 
output feedback controllers of flexible structures. Con- 
trollers for the Middeck Active Control Experiment have 
been based on LQG and Tii theory along with an ap- 
proach to maximum robustness to real uncertainty using 
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Popov multipliers [12]. Hoo and fi controllers are synthe- 
sized for the NASA - Langley Minimast Structure [4]. An 
approach to attenuate modal vibrations is researched us- 
ing eigenstructure assignment to control the eigenvectors 
that determine mode shape [26]. 

Full information feedback is an optimal feedback config- 
uration that provides direct measurements of all states 
and disturbances to the controller. Packard tt al demon- 
strate optimal full information controllers may be com- 
puted to minimize a fi upper bound for systems with 
complex uncertainty [21, 22]. These results are extended 
to compute optimal full information controllers that di- 
rectly account for additional information provided by 
purely real uncertainty in the synthesis algorithm [3]. 

The synthesis procedure involves the solution to a linear 
matrix inequality (LMI). LMFs have received a great 
deal of attention in recent systems and controls litera- 
ture [1, 7, 8]. LMI minimization problems may be solved 
using standard convex optimization techniques. The El- 
lipsoid Method and the Method of Centers are readily 
suited to LMI optimization [7]. 

This paper considers the LMI minimization for synthesis 
of optimal full information controllers of a high order 
flexible structure. There are two main results presented. 

• A globally optimal n controller is synthesized for 
a flexible structure. 

• Large dimension LMI optimization is solved on a 
Cray supercomputer. 

Vibration attenuating controllers are computed for a 
flexible structure. It is desired to attenuate the response 
at 11 natural frequencies between 0 and 70 (rad/stc). 
The controllers must be robust to parametric uncertainty 
in the modal damping coefficients and input multiplica- 
tive uncertainty. A globally optimal full information con- 
troller is compared to an output feedback controller com- 
puted using D-K iteration. Standard workstations are 
impratical for the optimal controller synthesis for the full 
order model due to the large computational cost. A Cray 
supercomputer is able to compute the optimal controller 
in considerably less time using a vectorized algorithm. 
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Robustness with Real/Complex Uncertainty 

Define x £ R n ' as the vector of states, z £ R n * as the 
vector of errors, y £ R n » as the vector of measurements, 
d £ R n- as the vector of disturbances and u £ R n ° as 
the vector of control inputs. The state-space description 
of a linear time-invariant plant can be represented as 
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where A £ R n ’ xn ',5i £ R n - xn *,B 2 £ R"- xnc ,Ci £ 
R n * xn, ,C 2 £ R n * xn * , and the E matrices of appropri- 
ate dimensions. 


Real parametric uncertainty is allowed to enter the prob- 
lem as scalar or repeated scalar blocks. Complex uncer- 
tainty enters the problem as scalar, repeated scalar or 
full blocks. 

The function p is defined as 

Pa(AO : dct(I — A/A) — 0} 

A 

The value of p depends on the block structure of A. Up- 
per and lower bounds for p have been derived which uti- 
lize two sets of structured scaling matrices. These scaling 
matrices are similar in structure to the uncertainty block 
structure. 


Define P as the Laplace transform of this system, 
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Define S(P) as the set of all real, rational, proper con- 
trollers, K(s ), which stabilize the closed-loop system. In- 
troducing a controller into the system leads to the fol- 
lowing linear fractional transformation (LFT). 

F,(P, K) = P n + PnK{I - P n K)~ l P n 


V = {D = ET = diag (D* . . .D*D? . . .D°d\ . ..d‘ p ) 

: D? G C R '* R ',Df G C c '* c \d { € C} 

The set of scalings Q affect only the real parametric un- 
certainty blocks. 


g = [G = G* = diag{G i . . .G m 0 . . .0) : G,- G C* x *-} 
Define an upper bound for real/complex p [10, 27]. 

fi(M) < inf a ((DA/Z? -1 + jG) (/ + G 2 )*) 

GeQ 


This LFT represents the closed-loop transfer function 
with the lower loop of P closed with the controller K. 
Analyzing performance using the infinity norm leads to 
the following minimization problem. 


inf 

KCS(P) 


\\F,(P,K)\\ C 


This is an Tioo optimal controller synthesis problem 
which has been solved using state-space equations [9, 11]. 


The structured singular value, p, can be used to deter- 
mine robustness of the closed-loop system to structured 
modeling uncertainty and the achievable performance 
level in the presence of uncertainty. The uncertainty is 
allowed to include both complex and real parametric un- 
certainty. 


The uncertainty description is structured with two types 
of blocks. The blocks are repeated scalar or full block 
matrices. Let integers m,n,p define the number of real 
scalar, complex scalar, and complex full blocks. Define 
integers Ri, . . Rm such that the i 1 * repeated scalar 
block of real parametric uncertainty is of dimension 
Ri x Ri. Define similar integers Ci, . . . ,C n to denote the 
dimension of the complex repeated scalar blocks. The 
structured uncertainty description A is assumed to be 
norm bounded and belonging to the following set. 

A = {diag ( S*I Rl . ..6 R I Rm 6?Ic, . . . / c . A, ...A,) 

: «/* € R ,Sf G C, A. G C e ’ xe >} 


The structured singular value provides a measure of ro- 
bustness in the presence of the defined structured uncer- 
tainty. The V and Q matrices are allowed to vary with 
frequency for linear, time-invariant (LTI) uncertainty. 
These scalings are restricted to be constant when com- 
puting robustness in the presence of linear time- varying 
(LTV) uncertainty [24]. 

The objective of control design is to maximize robust 
performance which corresponds to minimizing p in this 
framework. An approach to output feedback control de- 
sign which minimizes complex p is called D — K itera- 
tion. This technique is used to synthesize output feed- 
back controllers for the flexible structure example. D—K 
iteration tries to achieve the desired robust performance 
objectives by integrating “Woo control design with com- 
plex p analysis [2, 9, 19]. D—K iteration alternately 
minimizes the complex p upper bound with respect to 
K ox D holding the other variable constant. This tech- 
nique has been applied with great success to a variety 
of aerospace applications despite a lack of guarantee to 
reach the global optimum [2, 14, 23]. A more detailed 
discussion of D—K iteration can be found in Refer- 
ences [2, 20]. An approach to minimizing the p upper 
bound for the full information case is presented in the 
following section. 
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Optimal Full Information Controllers 

An algorithm to compute the optimally scaled, Woo full 
information (FI) controller fdr a system with linear, time- 
varying (LTV) uncertainty is presented by Packard et. 
a/ [21, 22]. This algorithm considers the synthesis prob- 
lem with constant V scalings included to allow infinitely 
fast time variation in the complex uncertainty. This sec- 
tion outlines this synthesis procedure and presents the 
main theorem given in Reference [21]. 

The open-loop system considered is a continuous-time 
plant. The controller is provided with direct measure- 
ment of the states and disturbances and is called the 
full information plant. Define Pji as the full information 
plant. 


The closed- loop Woo norm condition is shown to be 
equivalent to a maximum singular value condition in- 
volving the constant matrix form of the discretized full 
information plant. The following theorem demonstrates 
this equivalency [21]. 

Theorem 0.1 Given the full information plant, Pji(s), 
and the constant matrix discrete-time plant, Pji, along 
with the set V of scaling matrices, then the following are 
equivalent. 

1. There exists D £ V and stabilizing K such that 

\\D^F l (P fi (s) l K)D-i \\ o0 < l 

2. There exists D £ V and stabilizing K along with 
real X = X ^ > 0 such that with Z = diag(X , D ). 

a (z*F,(P/i,K)Z-l) < 1 

The resulting constant matrix problem has two scaling 
matrices, D and X . D is the original constant matrix that 
scales the plant inputs and outputs. X is the symmetric 
positive definite solution matrix in the Riccati formula- 
tion of the Woo problem. It scales the inputs and out- 
puts of the discrete-time state equation. The X matrix is 
included in the singular value test to restrict the eigen- 
values of the discrete-time state matrix to be less than 
one to guarantee stability. 

Now perform a change of variables. Denote the entries 
of the discrete-time plant as {R,U, V,T) and introduce 
Q to replace K ( I + TI \)~ 1 in the closed-loop LFT for 
notational convenience. 

Ft(Pji,K)= R+UK(I + TK)- l V = R + UQV 

A variant of Parrott’s theorem is used to formulate two 
maximum eigenvalue conditions equivalent to the con- 
stant matrix singular value condition introduced in The- 
orem 0. 1 . The resulting eigenvalue condition is presented 
in Lemma 0.2. 



Lemma 0.2 Given the discrete-time plant, P /it formu- 
late the closed-loop LFT as Ft(P/i, A) = R + UQV and 
define U± such that [U U±] ts invertible and U t U l = 0. 
TAen 

inf <F 

0€5(P) 

Z62 

if and only if there exists Z £ Z such that 

\(UZ(RZ~ l R T - Z~ X )U. l) < 0 

The proof follows directly from Parrott’s theorem [16, 
21]. There is an additional eigenvalue condition involving 
the matrix V and a perpendicular component Vx. defined 
similar to (7x- The V matrix in the LFT of the closed- 
loop system with the full information plant is square and 
invertible so V±, is null. The maximum eigenvalue condi- 
tion utilizing this variable is automatically satisfied. The 
remaining maximum eigenvalue condition is a linear ma- 
trix inequality (LMI) with the scaling matrix as the free 
parameter. 

The Woo synthesis problem for the continuous-time full 
information system is thus reduced to a single constant 
matrix maximum eigenvalue problem using elements of 
the discrete-time full information system. The following 
theorem summarizes the synthesis procedure for systems 
with infinitely fast varying LTV structured, complex un- 
certainty [21, 22]. 


(z*(R + UQV)Z-±) < 1 


Theorem 0.3 Given n, -state full information plant, 
Pji, and the set V, define the following. 

1. the augmented scaling matrices Z 

~_f\X °1 x € R n ' xn ', X = X T > 0 \ 

£ - U 0 D J • D € V / 

2. matrices R, U 

p_ f (I + A)(I-A)~ 1 V2a(I - A)~ l B\ 

[ y/2Ci(I — A) 1 Eu + Ci(I — A) Bi . 

r 

[ J ?12 + Cl(/ — A ) $2 . 

3. matrix Uj. such that [ U U ±] is orthogonal 

Then, there exists a stabilizing controller matrix K and 
a constant D £ P such that 

||d»F/ (Pji, k) a 

if and only if the following convex set ts nonempty. 

\Z 6 2: A m « f UT (RZ R‘ - o 2 Z) UA < 0} # {«} 


The optimal controller gains may be computed by scaling 
the continuous-time plant with the constant T) matrices 
and utilizing standard Woo state-space algorithms given 
in Reference [9]. 
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This synthesis algorithm assumes complex uncertainty 
for the system. These results are extended to directly 
account for real parametric uncertainty in the control 
process [3[. Accounting for real uncertainty may sig- 
nificantly increase the number of free parameters and 
consequently the required computation time for synthe- 
sis. The next section considers the cost of synthesizing 
optimal controllers with respect to complex uncertainty. 
The same arguments apply to the cost of synthesizing 
optimal controllers with respect to real uncertainty. 


Implementation of LMI Synthesis 


The previous section demonstrated optimal synthesis for 
full information feedback reduces to a single linear ma- 
trix inequality. LMTs may be solved using standard con- 
vex optimization methods. This section investigates im- 
plementation issues and computational cost associated 
with solving the optimal full information synthesis prob- 
lem. The convergence characteristics of algorithms are 
not discussed but rather the cost at each iteration of the 
minimization is dicussed. 

There are many convex optimization algorithms that 
may be used to solve linear matrix inequalities. The El- 
lipsoid Method and the Method of Centers have received 
a great deal of attention in recent control literature. The 
Method of Centers is an interior point algorithm that has 
been shown to be very efficient for a variety of systems 
and control problems [1, 7, 8, 13, 18]. Computational ex- 
periments indicate this method suffers from the type of 
constraints introduced in the full information synthesis 
formulation. This paper will concentrate on the Ellipsoid 
Method. 

The Ellipsoid Method utilizes a subgradient vector to 
compute an optimizing direction for the function to be 
minimized. The subgradient is defined as follows [25]. 

Definition 0.4 The function f : C n — ► C has the set of 
subgradients, S, at x. 

S = {g : f{y) > /(*)+ < g,y-x> V y) 


The subgradient vector obeys rules similar to the gradi- 
ent vector and is equivalent to the gradient at differen- 
tiable points of the function. The use of the subgradient 
allows the Ellipsoid algorithm to be applied to minimiza- 
tion problems involving nondifferentiable functions. 

The Ellipsoid algorithm is most easily visualized in a geo 
metric framework [5, 6, 7]. This algorithm operates on an 
n-dimensional space to minimize a function with n free 
parameters. An ellipsoid is defined in the n-dimensional 
space within which the optimal solution is assumed to 


lie. Each iteration of the algorithm computes a new el- 
lipsoid of decreasing volume. The center point of each 
ellipsoid is tested each iteration for optimality. A subgra- 
dient condition at the center point is used to eliminate a 
portion of the current ellipsoid that cannot contain the 
optimal point. A new ellipsoid is then computed such 
that its volume is the smallest possible while still con- 
taining the region of the old ellipsoid that may contain 
the minimizing point. The optimafsolution, when found, 
is the center of the last computed ellipsoid. 


The initial ellipsoid for the algorithm must, of course, be 
chosen such that it contains the optimal point. There is 
no known method to guarantee choosing an ellipsoid that 
contains this point. It has been shown in practice that 
selecting the initial ellipsoid to be a large ball centered 
at the origin is often a good choice. 

Many factors, including numerical conditioning of the 
matrices and choice of initial ellipsoid, greatly affect the 
convergence properties of the Ellipsoid Method. It is 
difficult to formulate an upper bound that is not overly 
conservative on the number of iterations required. Prob- 
lems of similar dimension have been shown in practice 
to require as few as several thousand or as many as sev- 
eral hunderd thousand iterations for convergence. The 
computational cost of each iteration is more tractable for 
analysis. 


Theorem 0.3 presented the LMI formulation for synthe- 
sis. This LMI involves minimizing the maximum singu- 
lar value of a matrix function by searching over a set of 
allowable scaling matrices. Consider the set, Z , of scal- 
ing matrices for a system with n, states and associated 
uncertainty scaling matrices V with dimension n u . 


Z = 


/ r x 

« ■ 

U 0 

D 


X € R n '* n ',X = X T > 0 
D ev e 


Every Z € Z is a block diagonal matrix whose blocks 
must be symmetric and positive definite. Symmetry is 
assured in the computation algorithm by only optimizing 
over the unique entries of each block. Restricting each 
block to be positive definite arises as a set of constraints 
of the LMI. 


The actual optimization problem to be solved is pre- 
sented below explicitly stating the constraints. Consider 
a system with n * uncertainty blocks and associated n* 
blocks in .the V scaling matrices. Define D\,D 2 , . . - , D n% 
as the blocks in the scaling matrices. 


min A (U[ 
zez v 

( RZR * - 

cx-Z) U L ) 

subject to 

Z„ = 

X 

> 0 


Z t = 

D 1 

> 0 


= 

D n . > 0 
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The total number of free parameters for this calculation 
can be calculated using the number of states, n Sy and 
the parameters for the dimensions of each type of block, 
Ci , c, and R*. Consider there are ±n(n + 1) free param- 
eters for a real, symmetric matrix of dimension n. The 
total number of free parameters is given as A . 

, r n s (n, + 1) , o Ri(Ri + 1) y- Ci(Ci + 1) 

N = 2 + ^ 2 + 2 +P 

The term involving the number of states, n,, is the num- 
ber of free parameters in the symmetric A scaling matrix. 
The first summation term computes the number of free 
parameters in the V matrices for the real parametric un- 
certainty blocks. The second summation term computes 
the number of free parameters in the V matrices for the 
complex parametric uncertainty blocks. The final term 
considers one free parameter for each scalar V matrix to 
commute with the complex full block uncertainties. 

The main cost of each iteration in the Ellipsoid Method 
is due to computing eigenvalues. The function to be min- 
imized is a maximum eigenvalue that must be computed 
at each iteration. Additionally, the maximum eigenvalue 
of each constraint matrix must be computed at each it- 
eration. The Ellipsoid Method computes a subgradient 
associated with the function that is the most violated, ei- 
ther the main function to be minimized or the constraint 
matrices. 


The eigenvector associated with the maximum eigenvalue 
must also be computed. A subgradient may be formu- 
lated utilizing this vector. The following lemma demon- 
strates the subgradient calculation for the function to be 
minimized. 

Lemma 0.5 Define the function /(A) = A (A m XA) for 
a given matrix A G C n * n and matrix X in the set of 
Hermitian matrices X — {A : X = X * G C nxn }. De- 
note v G C n as an eigenvector corresponding to the max- 
imum eigenvalue of f(X ). The following g is a suitable 
subgradient for /(A) for some unit vector rj G C n . 

g = AvrjT f v*A m 


Proof Use standard linear algebra properties to show the 
following for A, A G X and the inner product defined as 
< A,B>= Tr(B m A). 


\(A*(X + X)A) = 


> 

> 

> 

> 

> 


max u*(A* XA)u + u*(A* XA)u 

u*u = 1 

u m (A m XA)u + u m (A m XA)u 
T] m V* (A* X A)vrj + T) m V* (A* X A)vrj 
\(A'XA) + (r 1 *v m )(A*XA)(vT 1 ) 
\{A*XA) + Tr (X(Avw m v*A m )) 
A(A*AA)+ < < 7 , A > 


Define y = A *f A and x — A to show 


The matrix dimensions can be calculated for each eigen- 
value computation. Consider the dimensions of the R 
and U± matrices presented in Theorem 0.3 for scalars 
n , 7 the number of states, n C} the number of control sig- 
nals and n u , the number of uncertainty signals. 

R G + ) f/j_ £ j^(n f +n.)x(n f +n m -n c ) 

The dimension of the function whose maximum eigen- 
value is to be minimized can be computed. 

Ul (RZR* - a 2 Z) U L G R( n -+ n *-"‘> x ( n '+ n '‘- flc) 

The dimension of the matrix in the first constraint func- 
tion is the number of states, n Jt and the dimension of 
the matrix in the remaining constraint functions is sim- 
ply the dimension of each block in the scaling matrix, 
Ci ,Ci or Ri. 

Computing the maximum eigenvalue for the minimiza- 
tion function is generally going to be the largest com- 
putational cost at each iteration. It is easy to see this 
dimension can become quite large for plant with a large 
number of states and many uncertainties. The polyno- 
mial or exponential nature of most algorithms indicates 
the computational cost for this eigenvalue problem makes 
this optimal synthesis impractical for high order systems 
on standard workstations. 


A (y) > A(x)+ < g,y- x > 

which satisfies the definition of a subgradient. □ 

Most of the common operations performed in the al- 
gorithm are standard matrix-vector and matrix-matrix 
computations. Several software packages and libraries 
implement these as Level 2 and Level 3 BLAS subrou- 
tines. These routines may be highly optimized for a par- 
ticular architecture. The efficiency of Level 2 and Level 3 
BLAS routines can often be enhanced by implementing 
vectorized code. Unlike most workstations, Cray super- 
computers are designed for vectorized algorithms. The 
LAPACK package directly implements vectorized BLAS 
routines. This paper utilized the LINPACK, EISPACK 
and LAPACK packages for code generation. 

Analysis of the matrices utilized in the synthesis algo- 
rithms demonstrates there is no apparent structure or 
sparcity that may be exploited. Previous LMI optimiza- 
tions have been performed with thousands of free pa- 
rameters on a standard workstation by exploiting signif- 
icant structure in the computations. This algorithm will 
not demonstrate those qualities for a general engineer- 
ing problem. The problem dimensions can become large 
enough that this algorithm is impractical for a standard 
workstation. This paper utilizes a Cray supercomputer 
for solving a large order synthesis problem. 
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Flexible Structure : Modei 


linear working range of ±.5 inch. Control elements of 
Bay 3 are given in Figure 2 


Controllers are designed for vibration attenuation of an 
experimental flexible structure. The flexible structure in 
this example is constructed at the Dynamics and Con- 
trols Laboratory in the Department of Aerospace Engi- 
neering and Mechanics at the University of Minnesota. 
The structure is designed to place 12 lightly damped 
modes between 0 and 100 rad/ sec. The flexible structure 
is given in Figure 1. 



Figure 1: University of Minnesota Flexible Structure 

The structure consists of a rigidly held fixed top plate 
and four hanging bays. Each bay contains an aluminum 
plate. The top four plates are spaced .62 m apart while 
the bottom plate is .47 m below the third bay. Plates in 
Bay 2 and Bay 4 are triangular frames while the other 
plates are solid. Each bay is connected by thin hollow 
rods at the corners of the plates. All components of the 
structure are aluminum. 

Bay 3 contains sensors and actuators for control. Di- 
agonal rods connnect the top and bottom plates of this 
bay. Colocated linear force actuators and displacement 
sensors are placed along the diagonal rods. The force ac- 
tuators are voice-coil type actuators produced by North- 
ern Magnetics with a limit of ±2 pounds of force. The 
displacement sensors are Trans- Tek 0242 sensors with a 



Figure 2: University of Minnesota Flexible Structure - 
Bay 3 

Experimental transfer functions are formulated from the 
actuators to the displacements by commanding sinusoids 
of varying frequency to the actuators. System identifica- 
tion using curve fitting techniques and model reduction 
via balanced reduction computes a 22 th order model [17]. 
Transfer functions for the model from actuators to the 
displacement sensors are given in Figures 3,4 and 5. 





<r«»>«ncy <i«VWC) 

Figure 3: Transfer Function from Actuators to Displace- 
ment Sensor 1 for 22‘ A Order Model 
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Figure 4: Transfer Function from Actuators to Displace- 
ment Sensor 2 for 22** Order Model 



Figure 5: Transfer Function from Actuators to Displace- 
ment Sensor 3 for 22 th Order Model 
J 


Flexible Structure : Control Objectives 


It is desired to formulate controllers to attentuate vibra- 
tions and lower the peak gains of the open- loop system. 
The peak gains should be reduced approximately by a 
factor of 5 for modes less than 70 rad/sec. Performance 
weightings are included on the displacement measure- 
ments to account for the desired levels. 


Krj = 


1.2 


W. 


perj 


1 

1.2 


W 3 , 

rr perJ 


1 

2.0 


Input multiplicative uncertainty is included to account 
for unmodeled dynamics and neglected high frequency 
modes. A dynamic weighting, W mu it, is included on each 
actuator input channel. 

Wmu “ = l °(s+ 1000) 


This weighting function puts 1.0% uncertainty at low fre- 
quency and 100% uncertainty at 100 rad/sec with 1000% 
uncertainty at high frequencies above 1000 rad/sec. 


Parametric uncertainty is included in the plant model 
to account for inaccuracies in the modal damping co- 
efficients provided by the system identification algo- 
rithm [4]. Consider the state space representation of a 
second order system. 


1 

s 2 + s2uC + uj 2 


• 0 

1 

0 ■ 

2 

— U> 

-2u< 

1 

1 

0 

o J 


An additional input signal and output signal is included 
in the plant model to account for uncertainty in the 
damping coefficient £. A weighting W is included to scale 
the magnitude of the uncertainty parameter <5^ which is 
norm bounded by 1 in the f* framework. 



o 

1 

0 

0 ‘ 

1 

-u, 2 

-2wC 

i 

i 

s 2 + s2w<;(l + W6() + u> 2 

0 

1 

-'lutyr 

0 

0 

0 

IT 
0 _ 


Damping uncertainty is only included on the modes near 
60 rad/sec in the range of the performance weighing. 
These 4 modes have 50% uncertainty in the damping 
coefficient while the remaining low and high frequency 
modes have no damping uncertainty. The frequencies of 
the 11 modes and the range of damping coefficients are 
given in Table 1 
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Table 1: Modal Frequency and Damping Values 


The control signal is affected by a disturbance input. A 
constant weighting of Wdi%t = -5 is included to normal- 
ize the disturbance signal affecting each control channel. 
Sensor noise is also included in the system to affect the 
displacement sensor measurements. Constant weightings 
of W noise = .5 is included to normalize the noise affecting 
each sensor measurement. 

The magnitude of the control signal is included as a per- 
formance error to limit the amount of control actuation. 
A weighting of W ac t = -2 is included to normalize the 
control measurement for each actuator. 
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The open-loop plant model with uncertainty blocks and 
weightings is given in Figure 6. 



Figure 6: Open- Loop Flexible Structure Block Diagram 


Flexible Structure : Control Design 

The system given in Figure 6 contains two types of un- 
certainty. The parametric block representing uncertainty 
in the damping coefficients is given as A^. The damp- 
ing coefficients at each mode are independent of the other 
modes so this block is a scalar uncertainty block with four 
uncertainty parameters. The multiplicative uncertainty 
in the control signal, A is a full block uncertainty 
matrix with 3 inputs and 3 outputs. The uncertainty 
block structure is A. 

A = diag{6 \ , $ 2 , 63, 64 f A mti /g } 

The modal damping coefficient uncertainty parameters, 
are real parameter variations. Varying the 
magnitude of 6,* between ±1 implies the damping coef- 
ficient Q for the i ih mode varies between the maximum 
and minimum values given in Table 1. They are treated 
in the synthesis procedure are complex variations but 
they are considered as real parameters in the analysis. 
The multiplicative input uncertainty contains magnitude 
and phase information and is treated as a complex linear, 
time-invariant uncertainty. 

The uncertainty block structure is augmented with a per- 
formance block for computing robust performance. The 
performance block is a full block uncertainty matrix with 
6 inputs (3 disturbance and 3 noise signals) and 6 outputs 
(3 errors and 3 actuator penalties). This block is always 
treated as a complex linear, time-invariant uncertainty. 

An optimal full information controller is computed for 
this system. This controller, A>/, is synthesized as- 
suming all uncertainty parameters are complex, time- 
varying. A future project will provide computational 


resources for synthesizing a controller that directly ac- 
counts for the additional phase information provided by 
the real uncertainty parameters. 

An output feedback controller is also designed for this 
system using D- K iteration. The feedback measurements 
to the controller are noisy signals from the displacement 
sensors. This controller, denoted Kdk , is of order 34 
after a single D-K iteration. Kdk assumes all uncer- 
tainty is complex time-invariant. The Bode plot of Kdk 
is given in Figure 7. 




Figure 7: Bode plot of controller Kdk 

The closed-loop transfer functions for the controllers are 
given in Figure 8. Each controller attenuates the re- 
sponse at the natural frequencies. The full information 
controller is an optimal controller and achieves better at- 
tenuation levels than the output feedback controller. The 
low frequency performance of Kdk is especially poor. 



Figure 8: Closed- Loop Transfer Functions 
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Flexible Structure : Robustness Analysis 


Robust stability with respect to time-invariant uncer- 
tainty may be analyzed for the linearized plant model 
for each controller using the /i-Analysis and Synthesis 
Toolbox [2]. The upper and lower bounds for robust sta- 
bility fi with respect to real and complex uncertainty are 
given in Figure 9. The n values achieved are .90 for I<dk 
and .93 for Kpj. 



Frequency (r»<V*#c) 


Figure 9: Real/Complex Robust Stability /i 

The fi values for robust stability with respect to all com- 
plex uncertainty are similar to Figure 9. Treating the 
parametric damping uncertainty as complex does not de- 
crease the robustness of either controller. 

Robust stability for each controller is driven by the un- 
certainty at the higher modes. The peak fi values oc- 
cur at frequencies of the modes which contain damping 
uncertainty and a significant amount of multiplicative 
uncertainty. Robust stability is achieved at low frequen- 
cies where the uncertainty is low and also at frequencies 
higher than the natural frequencies. 

Kqk achieves a lower robust stability p value than the 
optimal controller Kfi- This does not contradict the 
synthesis theory. Kfi is designed to minimize robust 
performance and does not directly consider robust sta- 
bility. 

Nominal performance is also computed for both con- 
trollers. The p values achieved are 1.18 for I<dk an< ^ 
.46 for Kfi . The output feedback controller is unable 
to meet the performance specifications at low frequency 
while the full information controller achieves the desired 
performance goals over all frequencies. The nominal per- 
formance /i values agree with the poor low frequency 
performance shown in the closed-loop transfer function 
in Figure 8. The nominal performance p values are given 
in Figure 10. 



Fr*qu*ncy |r«y*«e) 


Figure 10: Nominal Performance p 

The robust performance p for each controller is given in 
Figure 11 for mixed real and complex uncertainty. The 
robustness levels are not noticeably changed by consid- 
ering the damping uncertainty as complex or real. The 
robust performance levels achieved are 1.23 for Kdk an ^ 
1.09 for Kfi • 



fntjjmncf (r»cW*c) 


Figure 11: Real/Complex Robust Performance fi 

Figure 11 shows the full information controller is able 
to achieve a robust performance /i value less than the 
output feedback controller, fj for Kfi is approximately 
20% less than Kdk with complex and real uncertainty. 

Neither controller is able to meet the desired robust- 
ness and performance specifications for the amount of 
included uncertainty. Kfi shows /i close to the desired 
levels and is only 9% from the desired goal. A reduction 
to 1/1.09 or 92% of the amount of input uncertainty re- 
sults in Kfi achieving robust performance. 

Robust performance for each controller is not affected by 
restricting the damping uncertainty values to be purely 
real. The output feedback controller, Kdk* * s driven 
by the additive noise which affects the achievable perfor- 
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mance at low frequency. The full information controller, 
AV/, is driven by the input multiplicative uncertainty 
at the high frequency modes. The additional damping 
uncertainty increases the level of difficulty in designing a 
controller but restricting that uncertainty to be real does 
not increase robustness for the controller. It is unclear 
that including real uncertainty in the synthesis algorithm 
would lead to more robust controllers. 

Time simulations of the open-loop and closed-loop sys- 
tems for each controller for the linearized plant model 
are given in Figure 12. The full information controller 
achieves the desired attenuation while the output feed- 
back controller shows slightly less attenuation. 



Figure 12: Time Simulations 


Flexible Structure : Computational Cost 

J 

The computational cost of computing the optimal full in- 
formation controller is dependent of the number of free 
parameters in the optimization. This number is deter- 
mined by the number of states and the uncertainty struc- 
ture in the set Z in Theorem 0.3. 

The open-loop system given in Figure 6 shows the block 
diagram for control synthesis. This system contains 
25 states from the flexible structure and weighting ma- 
trices. The uncertainty block structure contains 4 scalar 
blocks to represent the parametric uncertainty and a full 
block for the multiplicative uncertainty. The total num- 
ber of free parameters is N . 

N = ^25(25 + 1) + 4 + 1 = 330 

Each iteration of the optimization algorithm computes 
several eigenvalues and eigenvectors. The largest matrix 
used for the eigenspace computation is the matrix func- 
tion of dimension 38 given in Theorem 0.3. 

U' L (RZR* - a 2 Z) U L G R 38 * 38 


A single processor Cray-XMP supercomputer is used to 
compute the optimal controller in this paper. The al- 
gorithm implements a bisection search for the lowest 
achievable /i upper bound. The number of iterations 
required for each /j level are shown below. 


level 

iterations 

result 

1.10 

0.75 

0.92 

330261 

335701 

394868 

controller exists 
no controller exists 
no controller exists 


Standard workstations were initially utilized for con- 
troller synthesis since the eigenspace dimension is not 
extreme. The algorithm was terminated after several 
days without reaching an optimal controller. Each itera- 
tion with the eigenspace computation took only 1 second 
approximately but this implies 11 days are required to 
compute the million iterations needed for the controller 
synthesis. 

The Cray supercomputer is able to directly utilize vector- 
ized loops in the eigenspace computation to significantly 
reduce the execution time for each iteration of the op- 
timization. The execution time to compute the optimal 
controller was approximately 6 hours. 

Conclusion 

This paper considers a computational algorithm for op- 
timal controller synthesis. The formulation of a glob- 
ally optimal controller for full information feedback is 
shown to reduce to an LMI optimization. Implemen- 
tation of this optimization using the Ellipsoid Method 
requires many expensive eigenvalue and eigenvector com- 
putations. A flexible structure is presented as an exam- 
ple to synthesize controllers for vibration attenuation. 
The resulting problem formulation is of large enough di- 
mension with no apparent sparcity or structure that a 
standard workstation requires several days for the com- 
putation. A vectorizing Cray supercomputer is used for 
synthesis. 
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